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ABSTRACT 

The strong variability of blazars can be characterized by power spectral densities (PSDs) and Fourier 
frequency-dependent time lags. In previous work, we created a new theoretical formalism for describing 
the PSDs and time lags produced via a combination of stochastic particle injection and emission via 
the synchrotron, synchrotron self-Compton, and external Compton (EC) processes. This formalism 
used the Thomson cross section and simple J-function approximations to model the synchrotron and 
Compton emissivities. Here we expand upon this work, using the full Compton cross section and 
detailed and accurate emissivities. Our results indicate good agreement between the PSDs computed 
using the J-function approximations and those computed using the accurate expressions, provided 
the observed photons are produced primarily by electrons with energies exceeding the lower limit of 
the injected particle population. Breaks are found in the PSDs at frequencies corresponding to the 
cooling timescales of the electrons primarily responsible for the observed emission, and the associated 
time lags are related to the difference in electron cooling timescales between the two energy channels, 
as expected. If the electron cooling timescales can be determined from the observed time lags and/or 
the observed EC PSDs, then one could in principle use the method developed here to determine the 
energy of the external seed photon source for EC, which is an important unsolved problem in blazar 
physics. 

Subject headings: BL Lacertae objects: general — quasars: general — radiation mechanisms: non- 
thermal — galaxies: active — galaxies: jets 


1. INTRODUCTION 

Blazars, active galactic nuclei (AGN) with jets aligned 
with our line of sight, relativistically beam nonther- 
mal radiation towards the Earth. This emission ex¬ 
tends across the electromagnetic spectrum, from radio 
to 7 rays. Although the observed radio emission proba¬ 
bly comes from the superposition of many self-absorbed 
sync hrotron comp onents spread out over several par¬ 
secs (|Konig]iri98lll . the location along the jet where the 
higher-frequency emission originates is a matter of some 
controversy. This is particularly true for the most power¬ 
ful class of blazars, the flat spectrum radio quasars (FS- 
RQs) with strong broad emission lines. These objects 
are thought to emit y-rays primarily through Compton 
scattering of an external radiation field (known as ex¬ 
ternal Compton or EC). However, the primary external 
radiation field for EC is no t known, and could be pro- 
vided by the accretion disk (|Dermer fc_SchIickeisei] []^93l 
Hool), the broad line region (BLR; iSikora et al.lIlQQ^. 
or th e dust torus (jKataoka et al.lll999HBlazeiowski et al.l 
I 2 OOOII . This uncertainty in the seed photon source reflects 
the uncertainty in the location along the jet of the pri¬ 
mary 7 -ray emitting region, which one would expect to 
be within ~ 0.1 pc from the black hole (BH) if the scat¬ 
tering of disk radiation dominates, or between ^ 0.1 and 
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~ 1 pc if the scattering of BLR emission dominates, or 

1 pc if the scattering of dust torus emission do minates 
ijChisellini fc Tavecchioll2?)?)^ ISikora et al.ll2009ll . 

The association of 7 -ray flar es with the ejection o f su¬ 
perluminal components (e.g.. iMarscher et al.ll20l'^ in¬ 
dicates that the flares are coincident with the 43 GHz 
core, which is likely located at > a few pc from the black 
hole, outside the BLR (although see iNalewaiko et al.l 
l2014fl . The detection of > 100 GeV 7 rays from FS- 
RQs also indicates the 7 rays must originate from in¬ 
side the BLR to avoid 7 7 absorption with BLR photons 
(e.g.. lAleksic et alll2011[l . However, rapi d 7 -ray variabil- 
ity ob served in FSRQs s uch as 3C 454 . 3 (iTavecchio et al.l 
HoiS), PKS 1510-089 (|Browul [2flT3l : ISa.ito et aLI 1201 311 
and 4 C 21.35 (also known as PKS 122A21: IAleksic et al.l 
l2011f l limits the size of the emitting region, and, assum¬ 
ing the emitting region takes up the entire cross section of 
a conical jet, implies that it should be with in ^ 0.1 pc of 
the bl ack hole, and hence inside the BLR. iDotson et al.l 
lj2012li point out a way to distinguish between these sce¬ 
narios based on the different cooling timescales at two 
energies if the seed photons are from the BLR (e.g. with 
dimensionless energy for Lya cq = Eo/meC^ = 2 x 10“®) 
or from the dust torus (e.g. ^ 1000 K dust producing 
photons with dimensionless energy cq = 5 x 10“^) due to 
the variation of the Compton cross-section at high ener¬ 
gies. 

Blazars are bright and highly variable at all wave- 
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lengths. Their variability appears stochastic, charac¬ 
terized by power-law colored noise in their power spec¬ 
tral densi ties (PSPs') at all wavelengths, includin g ra¬ 
dio (e.g., iTrippe e t al.l 120111 Park fc Trippd 1201^. op- 
tical le.g.. iCh atteriee et al.l ^ 12t iWehrle et al.l 1201.11 


Edelson et al. 1201^ Revalski et alJ 1201411. X-rays (e.g. 


Zhang et al.lll999[ Kataoka et al.ll2001: Zham 


and 


7 rays (e.g., Aharonian et al.l l2007t lAbdo et all IMnI: 

iNakagawa &: Moril2013ll . Despite the popularity of PSDs 
for characterizing blazar variability, until recently a the¬ 
oretical motivation has been lacking. Previously, we 
presented a new theoretical formalism for the model¬ 
ing a nd interpretation of these sources (jFinke fc Becked 
12014 hereafter Paper I). This work was based on solv¬ 
ing the Fourier-transformed electron transport (continu¬ 
ity) equation for the case of synchrotron and Thomson 
cooling, yielding a solution for the electron distribution 
as a function of Fourier frequency. This solution was 
then combined with simple J-function approximations 
for synchrotron and Compton scattering emission to pre¬ 
dict the expected PSDs and Fourier frequency-dependent 
time lags one would expect to observe from blazars. We 
showed that these Fourier transform-related data prod¬ 
ucts contain detailed information about the characteris¬ 
tic timescales for particle escape, energy losses, and light 
crossing in the emission region, which is taken to be an 
outflowing spherical pla sma blob. Several o f thes e fea¬ 
tures were also noted bv IMastichiadis et al.l (1201311 . 

Here we extend the work presented in Paper 1, by 
solving the electron continuity equation with the full 
Compton cooling rates, including Klein-Nishina effects 
(Section [5]). We then compute the synchrotron and 
EC PSDs, going beyond the simple 5-function approx¬ 
imations adopted in Paper 1 and employing well-known 
precise expressions for the synchrotron emissivity and 
Compton scattering cross section (Section [3]). We also 
improve on our previous work by incorporating the light 
travel time effects specific to a spher i cal ge ometry, as de¬ 
scribed by iZacharias fc Schlickeiserl (|2013[1 . We do this 
in the context of EC emission from FSRQs. In Section 
|4]we compute the associated time lags, and in Section | 6 ] 
we describe how the electron cooling timescales, obtained 
from spectral breaks of PSDs, can be used to constrain 
the seed photon energy of the external radiation source 
for EC. This metho d is very close to the one outlined by 
iDotson et al.l ()2012f l. except that we suggest using PSDs 
to obtain the cooling timescales, rather than directly ob¬ 
serving the decaying part of the flare light curve. 

Note that this paper deals with the PSDs and time 
lags from nonthermal emission from relativistic jets. The 
PSDs from coronae in Galactic black holes and AGN 
exhib it breaks correlated with the mass of the black 
hole (|McHardv et al.l [200611 . This is true for the emis- 
sion from the coronae of AG N with jets as well (e.g., 
IChatteriee et al.l I2009L 120111) . There is some theoreti¬ 
cal indication that coronal emission, reflected off of an 
accretion disk, can create Pourier frequency-dependent 
time lags that may be relat ed to black hole mass 
(|EmmanouloDOulos et al.l [20141 ). However, in this pa¬ 
per, we consider only nonthermal jet emission, beamed 
away from the jet, so that no such features are expected. 
There is no expected relation between black hole mass 
and PSDs and time lags in our model. Although there 
may be a connection between jet and coronal emission. 


we do not consider such an effect here. 

2. ELECTRON DISTRIBUTION IN THE FOURIER 
FREQUENCY DOMAIN 

Consider a spherical homogeneous “blob” of Thomson- 
thin plasma of radius R containing nonthermal elec¬ 
trons and positrons (both of which are hereafter re¬ 
ferred to as electrons) at highly relativistic energies, and 
a tangled magnetic field of strength B. The electron 
distribution is assumed to be isotropic and the num¬ 
ber of electrons between Lorentz factors 7 and 7 -I- ^7 
as a function of time t is Ne{’^]t). Electrons may be 
injected into the blob at a rate < 5 ( 7 , t), lose or gain 
energy at a rate 7 ( 7 ,!), and escape from the blob 
on a timescale tesdjyt). In this case, the evolution 
of Nph;t) is described by (e.g., IMastichiadis fc KirkI 
1995t Chiaberge fc Ghisellinilll999t iLi fc Kusunosell200Cll : 


Bottcher fc Chianell2002HChen et al.ll201lL I2012D 




dt 


We will further assume that 7 , tesc, and R are indepen¬ 
dent of t. In this case, taking the Fourier transform of 
both sides of Equation o gives 

- 2-KifNe{lJ) + ^[7(7)^e(7,/)] + 

= Qh,f), (2) 

where the tilde refers to the Fourier transforrrQ, / is the 
Fourier frequency, ui = 27r/ is the angular Fourier fre¬ 
quency, and < 5 ( 7 ,/) is the Fourier transformed source 
term. It was shown in Paper I that if 7 < 0 and tesc is 
independent of 7 then the solution to Equation ([5]) is 


Neil,f) = 


17(7)1 


dl' Qil'J) 


X exp 


te 


— lU! 


di‘ 


17(7") I 


( 3 ) 


We make the ansatz that the Eourier-transformed source 
terms is 

= Qo(///o)-“/V''i?(7;7i,72)i?(/;/i,/2) (4) 

where 

(a;, a, ) ~ I Q otherwise ^ ^ 

is the boxcar function. We showed in Paper I that the 
normalization constant 


Qo = 


27rAt(Li„j) 


mec’^G^I^ + If- 2IrIo + 


( 6 ) 


where G, A, A, and /q are integrals defined in Appendix 
B of Paper I, and (Tinj) is the root mean squared power 
injected in electrons over a time interval At, so that At 
is the duration of the light curve. 

We will take into account two energy loss (i.e., cool¬ 
ing) processes for 7 ( 7 ): synchrotron and EC. We neglect 

^ See Paper I for the definitions of the Fourier transform and its 
inverse used here. 
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adiabatic losses, since we assume R is independent of t. 
The synchrotron cooling rate is given by 


where 


-isyil) 


dcCTT 

inieC^ 




UB = 


Stt 


( 7 ) 

( 8 ) 


is the magnetic energy density, c is the speed of light, 
me is the electron mass, and ut is the Thomson cross 
section. The Compton cooling rate, valid in the Thomson 
through Klein-Nishina regimes, for the case of scattering 
of an external, isotropic, monochromatic radiation field, 
is given by 


-lEcil) 


ScaTUg 

SrrieC^eg 


GBMsil^o) 


( 9 ) 


where ug and eo are the energy density and dimension¬ 
less photon energy, respectively, of the external radiation 
field measured in the frame of the host galaxy, and 


Gbms{x) 


8^ l + 5x 
3^(1 -b4a;)2 
ln(l -I- Ax) 


X 


3 3 1 

-f T~9 

X Ax^ 


Ax /'2 1 I \ 

1 -I- 4a; \ 3 2x 8a;^ ) 

ln[l -I- 4a;] ln[4a;] \ 

2a; X J 


5 1 

2 X 


1 

X 


E 


{l + Ax)-^ 



(jBottcher. Mause. fc Schlickeiseii [TO^ . A simple and 
fairly accurate approximation for the Compton cooling 
rate is given bv IModerski et all (I2005f l. and it is 


-i EC,Mil) 


dcCTT 

iirieC^ 


Ug-f'^ Mg{A^eg) 


( 11 ) 


the observer’s frame, except for the quantities related to 
the external radiation field, eg and ug which are in the 
frame of the host galaxy, so that e'g = Teg and Ug ~ F^ug. 

Below we describe the observed flux in the time 
domain {F{e,t)) and in the Fourier frequency domain 
(F(e, /)) for synchrotron and EC for the simple 6- 
function approximations and for the more precise expres¬ 
sions, and compare them with each other. This uses the 
result for Nf.il,f) described in Section^ The PSD can 
be found from the Fourier transformed flux by 

5(e,/) = |F(c,/)|2 = F(e,/)F*(c,/) (14) 

where the asterisk denotes the complex conjugate. In 
Paper I, we calculated the synchrotron emission includ¬ 
ing light travel time effects in a very simple, cylindri¬ 
cal geometry. Here, we include light travel time effects 
consistent with spherica l geom etry, using the results of 
iZacharias fc Schlickeiserl (|2013ll . See Appendix lAl for de¬ 
tails. 


3.1. Synchrotron 

3.1.1. Delta Function Approximation 

With the 5-function approximation, synchrotron pho¬ 
tons with observed energy e are created by electrons with 
a comoving Lorentz factor given by 


7 


/ 

sy 


e(l + ^) 

^Df-B 


(15) 


with corresponding total radiated power per electron (cf. 
Equation 0) 

* /■ \ 2 4cCT '(P 2 r -. 

Pe = -jsyiVmeC = —^UbJ ■ (16) 

The corresponding izFh synchrotron flux as a function of 
time t (including light travel time effects) if the blob has 
a tangled comoving magnetic field B is 


where 


Mg{x) 


1 

(1 + :e)3/2 ■ 


( 12 ) 


Naturally, the total cooling rate from synchrotron and 
EC losses is 


7 ( 7 ) = isyil) + iEcil) ■ (13) 


We neglect th e effects of SSC cooling, which can be quite 


complicated (ISchlickeised I200S 

; ISchiickeiser et al.l 120101: 

Zacharias & Schlickeiserl 120101 

2014 . 

2012ayH 120131: iZacha.ria.sl 


3. OBSERVED PSDS 

Now consider the case where the plasma blob described 
in Section [5] is moving with small angle 0 <C 1 to our 
line of sight at a speed that is a fraction j3 that of light 
(in the host galaxy frame) giving it a Lorentz factor F = 
(1-/3^)“^/^ and a Doppler factor 5 d = [r(l —/3cos0)]“^. 
The plasma blob comes from a supermassive black hole 
in a galaxy at a cosmological distance with cosmologi¬ 
cal redshift z. Primed quantities refer to the frame co¬ 
moving with the blob, while unprimed quantities refer to 




(SKsy,sil + z) 


higtic 


t'c 


t'c 


( , tSo 
X (r.r 


-t' 


(17) 


where Vg is the comoving electron distribution, 

2i?'(l -1-z) 


t/c - 


cSd 


(18) 


denotes the light-crossing time for the spherical blob, 




(19) 


cb = B/Bcn and B^r = 4.414 x 10^3 Q 
Similar to Paper I, the Fourier transform of this is 


47ri/(l -|- z)R' 


/ (1 + ^)/ 


5d 


X < exp 


c5d 


- 1 


( 20 ) 
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3 . 1 . 2 . Exact Synchrotron Emissivity 

The vFi, synchrotron flux as a function of time t in¬ 
cluding light travel time effects is given by 


F^He,t) = 


6K,y{l + z) 


Sotic 


dt' 


t'c 

m' 


t'c 

m' 




tSi) 

m 


- 1' ) Res 


f 2e(l -|- z) 

( 21 ) 


where 


K — 

J:\sy — 


•\/35De(l + z)e^B 


4,TThd\ 


and 


p'K pOO 

Rcs{x) = ^ d9 Sind / dy 

^ J0 Jx/ sin 9 


( 22 ) 


(23) 


is the function from iCrusius fc Schlickeised (|1986ll . 
The corresponding comoving power is again given 
by Equation (HU). Accu rate approximations to 

Rns(x ) are given by IZ irakas hvili fc Aharonianl 
(l2007f). iF inke. De rmer. fc Bottcherl (i20()^ . and 
iJoshi fc Bbttcheinioni) . Following the procedure 
outlined in Paper I, the Fourier transformed synchrotron 
flux using this precise expression is given by 


F^y{eJ) = 


QKsyjl + z) 
2TYiftic6ri 


exp 


47 rz/(l -I- z)R' 


X / dy'iV' 7 


(1 + ^)/ 
fc 


cfc 

Res 


- 1 


/ 2e(l -|- z) 

VSesfcy'^ 

(24) 


3 . 2 . External Compton 
3 . 2 . 1 . Improved Delta Function Approximation 

In Paper I we used a ^-function approximation 
for Compton scattering that was only valid in the 
Thomson regime. Here we make use of a very 
useful ^-approximati on for Compton scattering by 
iModerski et al.l (l2005f) that is fairly accurate (to within 
about 10%) in the Thomson (dfcy^j^^eo 1) through 

extreme Klein-Nishina (KN; 4fc7^(^eo 1) regimes. 

We assume the external radiation field is monochromatic 
with dimensionless energy cq that is isotropic in the frame 
of the host galaxy. Including light travel time effects, this 
yields for the EC flux as a function of time 


jpEC,5e^\ &KecA^ + z) 

— Ek. — 

X Afo(4fc7^(;^eo)Af2(4fc7^(^eo) 
fIB! Ic 

X / dt! 

Jo 


t'c 


1 

1 (XI 
_1 

\2R') _ 




where 


Kec,s = 


JL 

Sttc?! 


C(TtMo7bc 


(25) 


(26) 


Mi(x) = (y) = -^ 

7 


J^dy y Jc{x,y) 
J^dy Jc{x,y) 


and 


M2{x) 


dhi{x) 

d\n{xMi{x)) 


(27) 


(28) 


The function Mi{x) is computed using the “Jones 
formula” Compton scattering kernel for is o tropic 
electron and photon distributions (jJonesI 119681 
[Blumenthal fc Gouldlll97(ill . 


Jeix, y) = 2w Intc -I- (1 -|- 2w)(l — rc)-!- 
1 {xw)"^ 


X -- 


2 1-1- xw 


(1-w) , 


(29) 


in which 


w = 


y 

x(l - y) ■ 


(30) 


In Equation p5|) one can find ^'ec from e by solving the 
equation 


fc7^pMi(4fc7^geo) 
l + z 


(31) 


numerically for ^'ec- The function Mi{x) has the asymp¬ 
totes 


\ x/3 a; <C 1 Thomson Regime 
0.691 a; ^ 1 Extreme KN Regime ‘ ^ ' 


This implies 


Iec 


fcoTsceo < 1 
(o.ra'i)iD 4fc7^peo > 1 


and 

M2{x) 


1/2 X <C 1 Thomson Regime 
1 X 1 Extreme KN Regime 


(33) 


(34) 


Fourier transformation of the approximate EC flux 
given by Equation (1251) yields 


zEC,5 


(e./) = 


(!Kec,8{^ + ^) 


N, 


^TTitlcSjo 

47 rz/(l -b z)R' 


e IeC^ ■ 


( 1 +^)/ 


fc 


(35) 


exp 


cfc 


- 1 


X Alo(4fc7/;(;;eo)M2(4fc7^f;eo) . 


3 . 2 . 2 . Exact Compton Cross Section 

We follow iDermer et al.l (|2009f) to calculate the flux 
from EC emission resulting from the scattering of a 
monochromatic external radiation field that is isotropic 
in the frame of the black hole and host galaxy. We again 
include light travel time effects for a spherical geometry 
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following IZacharias fc Schlickeised (|2013ri . This gives 

ODtlc Jo 




t'c 

1 — 

(N 

7.^ 

m ~ 

\2R’) _ 




X Jc ( 45D7^eo, 


1 + z 
ep + z) 

j'Sb 

for the lyFi, EC flux. Here 

3 C(TTe^(l + uq 


Kec = 


Andl 


where 


(36) 

(37) 

(38) 

(39) 

Recall that 72 denotes the high-energy cutoff in the elec¬ 
tron injection spectrum (see Equation U). The Fourier- 
transformed EC flux is then given by 


+ ^) / 1 , /i , 
7min rtc 


2(5d 


eeo(l + z) 


and 


Tmax T 2 


F^^{eJ) = 


6KecO + z) 
‘2lTlftlcSY> 


exp 


47 ri/(l - 1 - z)R' 


cJd 


- 1 


7'^ fo 


X Jc (4(5D7^eo, 


e(l + ^) 

7'5d 


(40) 


3.3. Numerical Results 

As described in Paper I, for compact sources, the light 
travel time effect (Appendix will not be noticeable. 
In this section, we remove this effect to more easily dis¬ 
play the comparison between the S approximation and 
the full calculation. Without light travel time effects 
(valid for <C C(5d[/(1 + z)]~0i the syn chrotr on PSDs 
in the (5-function approximation (Section 13.1.11) and full 
calculation (Section I3.1.2L respectively, can be written 
as 


S^y’\e,f) = 


6Ksy,sO + z) 
Sd 


M I ' 

1 syi 


Sd 


and 


S^y{e,f) = 


f6K,y{l + z) 

I <5d 






(41) 


(42) 


A comparison of synchrotron PSDs computed the two 
different ways is shown in Figure [1} Agreement between 



Fig. 1.— The synchrotron PSD for the 5-function approxima¬ 
tion (das hed curves; Equation |41p. full calculation (solid curves; 
Equation |42|h and broken power-law fit (dotted curves; Equation). 
Parameters are q = 2, a = 1, = 10^ s, (Tinj) = erg s“^, 

At = 1 year, 71 = 10^, 72 = 10^, R' = 10^^ cm, B = 1 G, 
r = 5 d = 30, UQ = 10“^ erg cm“^, eq = 5 x 10“^, and z = 1. At 
this redshift with cosmology (h, Qm, r2A)=(0.7, 0.3, 0.7) the lumi¬ 
nosity distance dL = 2 x 10^® cm. The observed photon frequency 
{meC^€lh) or wavelength (/ic/[meC^e]) is shown. Dashed-dotted 
lines indicate / = (3tcool)~^ for each curve where tcool < ^esc- 


the ^-function approximation and the full calculation is 
good, except for the 120 GHz curve. For this frequency, 
photons come from electrons with 7 ^ < 71 . Clearly for 
this range the ^-function approximation is not accurate. 
With an observer-frame cooling timescale dehned by 


^cool(^) 


1 z d'y' 

7 17(7')! 


(43) 


the PSDs with tcooi < ^esc exhibit breaks at / = 
[3tcooi(e)]“^. As tcooi gets closer and closer to tesc, 
and eventually gets longer than teso this approximation 
breaks down. For synchrotron, 7 in Equation (j43ll is 
calculated from e using Equation (fT^ . The full curves 
“wash out” the many minima of the (5-function approxi¬ 
mation, so that the PSDs more closely resemble broken 
power-laws. The break has a magnitude of 2, i.e., the 
break is from S{e, f) oc /““ to S{e, f) oc /“(“+2). 

Without light travel time effects (i?(, <C C(5d[/(1 + 
z)]“^), the EC PSDs in the (5-function approximation 
(Section 13. 2. ID and full calculation (Section 13. 2. 2D . re¬ 
spectively, can be written as 




6Kec,sO + z) 

<5d 


2 


X Mo(4(5D7£ceo)Ar2(4(5D7^c'eo) 


X 


Ne 



0 + z)f \ 

Sd ) 


2 


(44) 
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Fig. 2.— The E C P SD for the 5-function approximation (d ashe d 
curves; Equation m, full calculation (solid curves; Equation Ha), 
and broken power-law fit (dotted curves). Parameters are the same 
as in Figure ^ The observed photon energy (meC^e) is shown. 
Dashed-dotted lines indicate / = (3tcool)~^ for each curve where 

^cool ^esc • 


and 




f6KEc{l + zy " 

V 

C/- h+^)/ 

L. 7'2 n 




(45) 


Figure [2] plots the EC PSDs computed the two different 
ways. Here for EC, 7 in Equation (1431) is found from e 
by solving Equation (j31ll numerically. Minima features 
in the ^-function approximation are again washed out in 
the full calculation, and agreement is good, except for 
the 1 MeV curve, where photons are generated primarily 
from electrons with 7 ' < 71 . As with the synchrotron 
PSDs, in the EC PSDs for tcooi < tesc, breaks are seen at 
/ = [3tcooi]~^, and again this approximation is less valid 
as tcooi approaches tesc- 


4. TIME LAGS 

Time lags as a function of Fourier frequency are given 
by 


Ar(ea, £ 6 , /) = ^ arctan 


Yi{ea,ebJ) ' 

YR{ea,ebJ)_ 


(46) 


(e.g.. Paper I) where yR(ea,eb,/) and Yi{ea,€b,f) are 
defined by 


FieaJ)F*iebJ) = YRiea,ebJ)+iYiiea.ebJ) . (47) 


Equations (H5)) and (ITTl) above are combined with Equa¬ 
tions (|M|l and (|4n|l to compute the time lags with the 
Moderski approximation and full calculation for the EC 
time lags. The results are shown in Eigure [31 Note 
that light travel time effects will play no part in time 
lags, since they are energy-independent. Our numer¬ 
ical results indicates that the time lags in the limits 
/ < [27rtcooi(eb)]“^ and / < [27r<cooi(eb)]“\ are 

Ar(eo, £6, /) ~ [tcooi(e&) - Aooi(£o)] /3 (48) 



Fig. 3.— EC flux time lags. Solid curv es show the lags computed 
with the full calculation (Equation Hq]), dashed curves sho w th e 
lags computed with the 5-function approximation (Equation 1351 V 

for the full calculation, where tcooK^) is the new definition 
from this paper. Equation (1331) . The Moderski 6 function 
approximation however is about 15% higher in all cases. 
This is also true for synchrotron time lags, which are not 
shown. Clearly the S function approximations are not as 
accurate for time lags as they are for PSDs. 

In this pape r, we use the standar d convention that 
£b < £a (e.g., iKroon fc Becked I2014D . so that positive 
lags refer to hard lags (i.e., the hard channel lags behind 
the soft channel), and negative lags to soft lags (i.e., the 
soft channel lags behind the hard channel). Thus, Pig- 
ure|3]shows our model only reproduces soft lags. Note 
that this is in contrast to Paper I, where we did not use 
this convention. I n X-ray observation s of BL Lac objects, 
both soft lag s le.g. lZhang et al.ir2002[) and hard lags (e.g., 
lZhang|l2002D have been observed. Our model can repro¬ 
duce the soft lags but not the hard lags, contrary to our 
discussion in Section 6.2 of Paper I. This is because the 
continuity equation treated here only includes losses (due 
to synchrotron and EC emission) and does not include 
any energization processes. Inclusion of acceleration ef¬ 
fects in our model may be able to explain the hard lags 
if the acceleration timesc ale is shorter than the cooling 
timescale, as suggested bv IZhan ^ (120021) . We will explore 
this possibility in future work (Lewis, Becker, & Finke 
2015, in preparation). 

5. COMPARISON WITH OBSERVATIONS 

In Section 6 of Paper I, we compared our model with 
several PSD observations. In this section, we update 
some of our conclusions, based on the improvements to 
our model. Results from Paper I not described below are 
unchanged. 


5.1. Oscillatory Features 

The use of the accurate expressions eliminate the 
oscillatory structures in the PSDs seen with the 6- 
approximations, so that the PSDs more closely resem¬ 
ble broken power-laws. Possible minima from observed 
PSD, suc h as the VHE PSD o f PKS 2155-304 pre¬ 
sented b ylAharo nian et ^ (I2007D, or the QPO s obser ved 
by e.g. Lachowicz et al.l ( 20091) . IGupta et al.l (|2009D or 
IRani et al.l ( 2010i) 7 we now believe are unlikely to be as¬ 
sociated with the minima features predicted in Paper 1. 
These features seem to be an artifact of the (5-function 
approximation and they are not present in more accurate 
calculations. 
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5.2. The j-ray PSD of 3C 454-3 

As reported by INakagawa fc Morll (1201311 , the LAT 7 - 
ray PSD of 3C 454.3 has a break at a frequen cy of 1.5 x 
10“® Hz (although see lSobolewska et al.l201^ . If at LAT 
energies tcooi < ^esc, this implies a cooling timescale of 
2.3 X 10^ s in the observer frame. For parameters F = 
(5 d = 30, H = 1.0 G, uo = 10“"^ erg cm“^, 72 = 10®, and 
assuming the LAT is dominated by photons at 100 MeV, 
one can numerically solve for the seed photon energy, and 
obtain eg = 1-8 x 10 “®, which is quite close to the value 
for Lyof (2.0 x 10“® in units). For the same set of 
parameters, except Uq = 10 “^ erg cm“®, one gets eg = 
7.0 X 10“^ (or 0.36 keV), and for uq = 10~® erg cm“®, 
one gets eg = 3.3 x 10“® (or 1.7 keV). The latter two 
are unrealistically high for seed photon energies, and Lya 
does make more sense as a seed photon source, since it 
is the most prom inent line seen in optical spectra (e.g., 
iTelfer et al.ll200^ . however this is an unrealistically low 
energy density. In this regard, the Klein-Nishina effects 
do not alter our conclusions from Paper 1. 


5.3. The Optical PSDs of blazars 

lEdelson et al.l (I2013D have used Kepler observations 
to identify a “bend” in the PSD of the BL Lac 
W2R1926+42 corresponding to a period of 4 hours. As¬ 
suming the emission detected by Kepler is dominated 
by detections at 5000 A, and for similar parameters 
as in Section 15.21 with uq = 10“^ erg cm“® one gets 
eg = 2.3 X 10“^, about what one would expect for scat¬ 
tering of dust torus emission; for ug = 10 “^ erg cm“® 
one gets eg = 2.5 x 10“®, about what one expects for 
scattering Lya photons; and for ug = 10“^ erg cm“® one 
gets eg = 1.2 X 10“'*. The latter is probably an unre¬ 
alistic value for eg, and the values of ug for the former 
are about what one would expect for these seed pho¬ 
ton sources, although for Lya the energy density is a 
bit low. Thus for this source, scattering of dust torus 
emission seems the most likely, although this conclusion 
is strongly dependent on the assumed parameters. In 
Section [ 6 ] we describe a method of determining the seed 
photon source that is more model-independent. 

6 . METHOD FOR DETERMINING SEED PHOTON ENERGY 

Here we outline a technique for determining the energy 
of the seed photon source for EC. We assume the external 
seed photon source can be approximated as isotropic in 
the galaxy’s frame, and that it is monochromatic with 
dimensionless energy eg . We also assume from the breaks 
or time lags in the PSDs of blazars at several energy 
ranges in EC one can obtain the cooling timescale, as 
described in Sections [3] and SI For example, one might 
find breaks in Fermi-LAT PSDs at meC^Ca = 0.1 GeV, 
nieC^eb = 1.0 GeV, and mgC^ec = 10.0 GeV, and thus 
the cooling timescales at these energies. 

If one is observing synchrotron or EC, one can use 
Equation (S31) to compute the observer frame cool¬ 
ing timescale from synchro-Compton losses, which we 
rewrite as 


Aooi(c) — 


3(1 -I- z)mec'^eg 


dx 


1 


CU'tUb 


x'^ 1+ AcMg{x) 


(49) 


where 


and 


X2 = 4F72eg 



ub 


(50) 

(51) 


is the Compton dominance, assuming Jd = F. If the 
cooling timescale is estimated from a PSD or time lag 
that is emitting synchrotron, the integral’s lower limit is 


Xl,. 


= 4eg 


e(l -I- 2 :)fc 


es 


(52) 


If the cooling timescale is estimated from EC, then the 
integral’s lower limit is 

, 11 ^ / 4eeg(l -I- z) 

= M,g,.Ec] 

_ f 2y'3eeg(l -b z) dFy^^'eg < 1 /.on 

\ 4eeg(l + z)/(0.691) AT^'^ceg » 1 ' ^ ^ 


For computing xi above we have made use of the S 
function approximations for synchrotron and EC, as de¬ 
scribed in Section 13.1.11 and 13.2.11 For EC, a function 
created from three cooling timescales. 


r(eo,e&,ec) = 


Aool(Ca) ^cool(Cc) 
^cool (Ca ) Aool(cb) 


(54) 


is dependent only on eg and Aq- In principle, Ac can 
be determined from the broadband spectral e nergy dis¬ 
tribution fe.g.. lMever et al.ll2012l : lFinkell2013H . For syn¬ 
chrotron, this function is additionally dependent on the 
ratio Sd/cb- Observations of FSRQ synchrotron PSDs 
and lags will also suffer from contamination by thermal 
emission, from the accretion disk and dust torus. There¬ 
fore we will only be concerned with the 7 -ray emission 
from EC. 

We plot the r(0.1 GeV, 1.0 GeV, 10.0 GeV) for different 
values of Ac in Figure IH These plots demonstrate that 
given Ac and tcooi for three energies, one can estimate 
the energy of the seed photon source, eg. For low values 
of eo (eo ^ 10”®)5 the scattering at all three energies will 
be in the Thomson regime, and the value will be mostly 
independent of eg or Ac- Only at higher values of eg will 
Klein-Nishina effects become apparent. 

Although in principle it may be possible to deter¬ 
mine eg from PSDs or time lags computed from Fermi- 
LAT light curves, in practice this will be extremely 
difficult. During the brightest flares , it is possible 
to p r obe timescales ^ 1 hour (e.g . . iTavecchio et al.l 
20101: lAbdo et al.l 120111: iBrownI [20131 ISaito et al.l 120131 
Nalewaikoll2013ll . but this is rare. More typically, inte¬ 
grations over a few days to a week are need ed to signif¬ 
icant ly detect a source with the LAT (e.g., lAbdo et"^ 
IMS), so that probing the ^ 1 hour timescales necessary 
to detect the breaks is unlikely in most cases. So far, 
three FSRQs have published spectra from imaging at¬ 
mospheric Cherenkov tele s copes (lACTs) (|Albert_et_aL 
200 a lAleksic et all IWl |H ■E.S.S. Collaboration et ^ 
20W ), and one more has been detected by lACTs 
([Mirzov^ I 2 OI 5 I iMukherie^ 1201511 . It is certainly rea¬ 
sonable to expect t hat the proposed Cherenkov Tele¬ 
scope Array (CTA; lActis et al.l 1201111 will be able to 
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Fig. 4. — The function r(0.1 GeV, 1.0 GeV, 10.0 GeV) from Equa- 
tion ll54t plotted as a function of dimensionless seed photon energy 
eo for different values of Compton dominance Ac- 

get significant detections with shorter integration times 
than the LAT, and therefore it may be able to probe 
these timescales. However, at higher energies, the cool¬ 
ing timescales will also be shorter, making the breaks 
more difficult to observe. If the breaks are on sub-hour 
timescales, a time baseline (i.e.. At in our notation) 
of about 10 times the timescale of the expected break 
should be sufficient, so that a single night of observing 
would probably be a long enough baseline to observe 
a break. Probing short enough timescales would be a 
greater issue. Further, any source detected at > 20 GeV 
is almost certainly making 7 rays outside the BLR, oth¬ 
erwise the 7 -rays would be att enuated by 77 absorptio n 
with the Lya photons (e.g., iStern fc Poutaneiil 120141 1. 
However, an analysis with CTA could verify our results. 

The technique outlined here is similar to obtaining the 
cooling timescales by modelin g the decay t i mes o f in¬ 
dividual flares, as discussed bv IPotson et al.l (j 2012 tl . al¬ 
though our model has a few advantages. For example, 
with the new method, one is able to take advantage of a 
long series of data, rather than individual flares, possi¬ 
bly leading to more statistically significant results. When 
modeling individual flares, it is not clear if the decrease 
in flux from a flare is due to radiative cooling or due 
to the rapid decrease in the electron acceleration rate. 
The disadvantage of our technique is it may be that dif¬ 
ferent flares from the same object have different origins. 
Perhaps some flares occur close to the BH, so that the 


seed photon source is from the BLR, while others occur 
farther away, so that the seed photon source is the dust 
torus. PSDs based on long timescale light curves would 
have a difficult time disentangling these cases. 

7. DISCUSSION 

We have expanded upon our previous theoretical sce¬ 
nario for computing PSDs and time lags (Paper I). We 
have included the full Compton cross section in our solu¬ 
tion to the Fourier-transformed continuity equation. As 
in Paper I, we find breaks in the resulting PSDs asso¬ 
ciated with the cooling timescale, however, we must re¬ 
vise our definition of the cooling timescale. We have 
compared the simple (5-function approximations for syn¬ 
chrotron and Compton emission to more accurate ex¬ 
pressions. For PSDs the agreement between the approx¬ 
imations and the accurate expressions is very good, ex¬ 
cept for photons produced primarily by electrons with 
7 ' < 7 (. As in Paper I, we find that the breaks in syn¬ 
chrotron and EC PSDs always have magnitudes of 2, i.e., 
the break is from S{e, /) oc to S{e, f) oc 

At low Fourier frequencies, time lags calculated with 
the ^-approximations give lags about 15% greater than 
lags computed with the accurate expressions. The cause 
of the deviation is not clear. However, as in Paper I, we 
confirm that the time lags at low Fourier frequencies can 
be associated with the difference in cooling timescales of 
the two energy channels. 

Based on our theoretical work, it seems that one could 
in principle measure the cooling timescale based on ob¬ 
servations of blazar PSDs or time lags. If one can mea¬ 
sure the cooling timescale from 7 -rays produced by EC 
at three different energies, while also determining the 
Compton dominance from the blazar’s broadband SED, 
then in principle determine the energy of the seed pho¬ 
ton source. This would give a strong indication as to the 
location of the 7 -ray emitting region, something which 
has so far eluded understanding. Such a feat is unlikely 
to be possible with currently operating 7 -ray telescopes 
such as the Fermi-LAT, but may be possible with CTA. 

We are grateful to the referee for helpful comments 
that have improved this manuscript, and to C. Dermer 
and S. Larsson for useful discussions. JDF is supported 
by the Chief of Naval Research. 


APPENDIX 

A. LIGHT TRAVEL TIME DELAY 

In this section we derive the light travel time effect, directly following iZacharias fc Schlickeisej (I2013L c.f. Appendix 
C of Paper I). The spherical emitting region with radius R, volume V = AttR^/3 is homogeneous, with variations 
taking place throughout the sphere simultaneously. Its total emitting flux is F at time t, but the observer will see 
light from the closer part of the sphere arrive before the light from the farther part of the sphere. For the calculation 
of this effect, the emitting region is divided into an infinite number of “slices” each with infinitesimally small thickness 
dx emitting flux dF at time t. The slices are cut perpendicular to the direction of the observer so that each one has 
radius 


y = v^i ?2 _ (i^_a ;)2 ^ 

cross sectional area 


(Al) 


A = 7ry^ = tt{2Rx — x^) 


(A2) 
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Fig. 5.— Geometric sketch of spherical emitting blob for the purpose of computing light travel time effects. This sketch is in the frame 
co-moving with the blob. 


and volume 


dV = A dx 


(A3) 


(see Figure [5]). The flux emitted by each slice as a fraction of the whole is proportional to the volume of the slice as a 
fraction of the whole, i.e., 


dF _ dV 
~F ~~V ' 

The observer sees flux from each slice delayed by time t = x/c so that the observed flux at time tabs is 


(A4) 


Fobsitobs) = j dF{tobs -t)= J P^tobs - t) . 


3c 

T 




dt F{tobs - t) 



(A5) 


Combining this with standard formulae for synchrotron and Compton emission for F{t) fe.g.. lDermer fc Menonl[2009fl 
leads to the results in Sections lO and o 
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